Three-body problem in heteronuclear mixtures with resonant interspecies interaction 
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We use the zero-range approximation to study a system of two identical bosons interacting res- 
onantly with a third particle. The method is derived from effective field theory. It reduces the 
three-body problem to an integral equation which we then solve numerically. We also develop an 
alternative approach which gives analytic solutions of the integral equation in coordinate repre- 
sentation in the limit of vanishing total energy. The atom-dimer scattering length, the rates of 
atom-dimer relaxation, and three-body recombination to shallow and to deep molecular states are 
calculated either analytically or numerically with a well-controlled accuracy for various energies as 
functions of the mass ratio, scattering length, and three-body parameter. We discuss in detail the 
relative positions of the recombination loss peaks, which in the universal limit depend only on the 
mass ratio. Our results have implications for ongoing and future experiments on Bose-Bose and 
Bose-Fermi atomic mixtures. 

PACS numbers: 34.50.-s, 67.85. Pq 



I. INTRODUCTION 

The universal properties of particles with resonant 
short-range interactions are a subject of intense research. 
Such systems are characterized by a large scattering 
length and display universal phenomena associated with 
a discrete scaling symmetry [HE]- For identical bosons, 
Efimov found that there are infinitely many trimer states 
with an accumulation point at the scattering threshold 
when the s-wave scattering length a is tuned to the uni- 
tary limit 1/a — [5] : 

where m is the mass of the particles, sq ~ 1.00624, and 
IS the binding wavenumber of the branch of Efimov 
states labeled by n*. The geometric spectrum in ([T]) is 
the signature of a discrete scaling symmetry with scaling 
factor e^/'^° « 22.7. For a finite scattering length larger 
than the range of the interaction, the universal properties 
persist but there is only a finite number of Efimov states. 

Ultracold atoms are an ideal tool to study such phe- 
nomena since the scattering length can be tuned exper- 
imentally using Feshbach resonances. Efimov trimers in 
ultracold atomic gases can be observed via their signa- 
ture in three-body recombination rates 0HH]- Kraemer 
et al. provided the first evidence for Efimov trimers in an 
ultracold gas of ^'^^Cs atoms by observing the resonant 
enhancement of three-body recombination caused by the 
trimers [5]. In a subsequent experiment with a mixture 
of ^'^■^Cs atoms and dimers, Knoop et al. observed a res- 
onance in the loss of atoms and dimers [TU] which can be 
explained by an Efimov trimer crossing the atom-dimer 
threshold |llj . 

Several recent experiments have also obtained evidence 
of Efimov physics with other bosonic atoms. 



Zaccanti et al. measured the three-body recombination 
rate and the atom-dimer loss rate in an ultracold gas of 
•^^K atoms [T2]. They observed loss features at large pos- 
itive and negative values of the scattering length, posi- 
tions of which agree with the discrete scaling symmetry. 
Gross et al. measured the three-body recombination rate 
in an ultracold system of '^Li atoms [13] . They observed 
a three-atom loss resonance and a three-body recombina- 
tion minimum in the same universal region on different 
sides of a Feshbach resonance. Their positions are consis- 
tent with the universal predictions with discrete scaling 
factor of 22.7. Using ultracold ^Li atoms as well. Pollack 
et al. [H] observed 11 three- and four-body loss features 
in the inelastic loss spectrum. Their relative locations 
on either side agree well with the universal theory, while 
a systematic deviation from universality appears when 
comparing features across the resonance. The origin of 
this deviation is not understood. Barontini et al. |15| 
investigated the Bose-Bose mixture ^^K-^^Rb and found 
three resonance positions in the three-body loss. The two 
features for negative scattering length were attributed to 
the two possible Efimov trimers of the system, K-Rb-Rb 
and K-K-Rb, hitting the three-atom threshold. 

The Efimov effect can also occur for fermionic atoms 
with at least three spin states. The first experimental 
studies of many-body systems of ^Li atoms in the three 
lowest hyperfine states were recently carried out by Ot- 
tenstein et al. [IBj and Huckans et al. [T7]. Theoretical 
calculations of the three-body recombination rate sup- 
ported the interpretation that the narrow loss feature 
arises from an Efimov trimer crossing the three-atom 
threshold [T5ti2Ij . Very recently, another narrow loss 
feature was discovered in the region of much higher mag- 
netic fields by Williams et al. [22] and by Jochim and 
co-workers [l^. In this region, the scattering length is 
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much larger and several recombination features have been 
predicted using the universal theory |24j . 

In this paper, we focus on heteronuclear systems with 
two species of atoms where only the interspecies scatter- 
ing length is large. For comparable masses the scaling 



factor is quite large (for equal masses 



1986.1) 



as we now have only two resonant interactions out of 
three. However, in the case of two heavy atoms and one 
light atom, this factor can become significantly smaller 
than the value 22.7 for identical bosons [HIISHIT], which 
should stimulate experimental investigation of the dis- 
crete scaling invariance. Relaxation and recombination 
losses near an interspecies resonance have recently been 
investigated in mixtures of rubidium and potassium. 
The Bose-Fermi combination ^^Rb-'*°K has been stud- 
ied at the Joint Institute for Laboratory Astrophysics 
(JILA) [55] and measurements on the Bose-Bose mixture 
4ij^_87j^j^ have been carried out in Florence [TS]. We ap- 
ply our theory to these and other mixtures of interest for 
ongoing and planned experiments. 




FIG. 1: Integral equation for the atom-dimer scattering am- 
plitude A. Solid (dashed) lines denote atom species 2 (1). 
Mixed double lines denote the full dimer propagator. 



where P ~ \P\, fi — mim2/{mi + 771,2) is the reduced 
mass, M — rrii + 7712 is the mass of the dimer, and the 
limit of e — > +0 is understood. The dimer wave function 
renormalization is given by — g|a/i^/(27r). 

The scattering between a dimer and an atom is de- 
scribed by the integral equation shown in Fig. [T] Using 
the Feynman rules derived from Eq. ^ and projecting 
on relative s-waves, we have 



II. METHOD 

In this section, we set up the effective field theory 
method which provides a convenient implementation of 
the universal theory for large scattering length. We set 
H = 1 but restore the dimensions in our expressions for 
the recombination rate constants. We consider a system 
of one boson or fermion of mass mi (species 1) and two 
identical bosons of mass 7712 (species 2). We assume the 
interspecies interaction to be resonant and characterized 
by the s-wave scattering length a ^ £vdw, where £vdw 
is the van der Waals range of the potential. Nonreso- 
nant intraspecies interaction will be neglected. If species 
1 is also bosonic and weakly interacting, all the forth- 
coming results directly apply to the other possible (in- 
teracting) triple by simply exchanging the labels 1 and 
2. We therefore include only the interaction between the 
atoms of species 2 and the dimers. Hence, our effective 
Lagrangian reads 



V'i idt + 



ipi+tpl{ idt + 



V'2 



27771 / V 27712 



g2d^d 
, (2) 



where the dots represent higher-order derivative interac- 
tions, and (72 and 53 are the bare two- and three-body 
coupling constants. 

From the Lagrangian ([2]), we can deduce Feynman 
rules and obtain the full dimer propagator and the three- 
body integral equation (see Ref. P4 for details on the 
derivation). For the full dimer propagator we get 
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where (j-ad = 7772(7771 -I- 7772)7(27772 -t- 7771) is the reduced 
mass of an atom and a dimer, the relative momenta of 
the incoming and outgoing atom-dimer pair are denoted 
by p and A:, respectively; and E is the total energy. The 
contribution of the s-wave projected one-atom exchange 
is given by 



Kip,q) 



1 
2m 



In 



p^ + g^ + 2pg^-2^g-7£ 
p^ + q^ - 2pq:^^ - 2^iE ~ le 



(5) 



and the contribution of the three-body coupling 53 can 
be written as 



93 



g(A) 
A2 



(6) 



where H(A) is a dimensionless log-periodic function of 
the cutoff A, which depends on a three-body parameter 
A* [29]. The mass-ratio dependence of the discrete scal- 
ing factor exp(7r/so) follows from the equation for sq: 



So cosh(7rso/2) — 2 s'mh{(f)So) / sin(20) = 0, 



(7) 



M^"-2m)- 



(3) 



where we introduce the parameter 

(j) = arcsin [1/(1 + S)] 



(8) 



and the notation S = toi/to2- For particles of equal 
mass, the solution of Eq. ([t]) is sq ~ 0.4137 leading to 
the scaHng factor exp(7r/so) ~ 1986.1. Because of the 
log-periodicity of H{A) one can always find a value of 
the cutoff A with H — 0. In practice, one can therefore 
simply omit the three-body coupling in the leading-order 
calculations and use the cutoff A as a three-body param- 
eter [3^ . We use this strategy in the following. For fixed 
5, the values of A and A* are related by a multiplicative 
constant. 

The scattering amplitude A has simple poles at the 
three-body bound-state energies E = —Bt < 0. The en- 
ergies can be obtained from the solution of the following 
homogeneous integral equation for the bound-state am- 
plitude B: 

which has nontrivial solutions only for three-body bind- 
ing energies Bt > 0. In the following, we use Eqs. Q 
andQ to describe three-body properties of heteronuclear 
mixtures. 



III. NUMERICAL RESULTS 

Few-body loss phenomena offer a unique view on scat- 
tering processes in ultracold quantum gases. In particu- 
lar, an enhancement of the loss rate can be an evidence of 
a few-body resonance. The universal theory predicts the 
relative positions of such resonances as a function of the 
scattering length. The universality can thus be tested ex- 
perimentally by measuring the lifetime of a cold atomic 
gas as a function of a. Ideally, in order to see the univer- 
sal scaling, one needs to detect more than one resonance 
in a single universal region, that is, a region where the 
three-body parameter can be assumed constant. This is 
believed to happen in a narrow vicinity of a Feshbach 
resonance, where large variations of a are accompanied 
by (assumed) much weaker variations of the three-body 
parameter. We now discuss three-body loss resonances 
in a heteronuclear mixture as predicted by the universal 
theory. 

A. Resonance positions 

The mechanism of three-body losses and its relation 
to the positions of Efimov levels in the heteronuclear 
case are qualitatively the same as for three identical 
bosons. The scattering-length dependence of the en- 
ergy of a generic trimer is illustrated in Fig. |2] On the 
negative side of a Feshbach resonance the trimer hits 
the three-body scattering threshold at a = a_ < 0, 
which leads to an enhanced probability of finding three 
atoms at distances of the order of \a\. Such atoms can 
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FIG. 2: (Color online) Dependence of the trimer energy on the 
inverse scattering length 1/a in arbitrary units (dashed line). 
The parameters a_ and a, specify where the trimer state hits 
the three-atom and atom-dimer thresholds, respectively. 



then approach each other to distances of the order of 
£ydw s-nd recombine into a deeply bound dimer and a 
residual atom. The released binding energy [of order 
fi^ / {2iJ,£lj^^)] transforms into the kinetic energy of the 
recombination products, which hence leave the trap. On 
the positive side of the Feshbach resonance there exists 
a weakly bound (shallow) dimer state with binding en- 
ergy Bd = h? / [2^0?). This formula taken with minus 
sign determines the atom-dimer threshold (solid line in 
Fig. [2| . By following the dashed line in Fig. [2] from neg- 
ative to positive values of a, one can see that the trimer 
crosses the atom-dimer threshold at a = a* > 0, where 
one predicts an elastic atom-dimer resonance. At this 
point formation of deep dimer states (in this case called 
relaxation) in atom-dimer collisions is also enhanced for 
the same reason as above. According to [12], the atom- 
dimer scattering resonance should be noticeable even in 
a purely atomic sample due to rescattering processes. In- 
deed, before leaving the trap, shallow dimers formed in 
the process of three-body recombination can collide with 
other atoms. The recombination rate itself is featureless 
around a = a*, but the atom-dimer cross section in the 
vicinity of this point is highly a dependent. Thus, at 
a = a* the three-body recombination can be enhanced in 
the sense that many more than three atoms are expelled 
from the trap leading to a measurable trap loss. We come 
back to this issue in Sec. IVl 

The ratio of the two resonance positions, a*/|a_|, is 
of fundamental importance for studies of the universal 
three-body physics as in the universal limit it does not 
depend on the three-body parameter. In order to calcu- 
late this ratio we solve the bound-state Eq. ^ for Bt — 0, 
a < 0, and for Bt = Bd, a > 0, with the same (arbitrary) 

cutoff A. The solid line in Fig. [3] shows al"Vk^^l ^ ^ 
function of the mass ratio S. Here we use the index n in- 
troduced in Eq. ([T]) in order to emphasize that the values 
of a* and a_ are taken for one and the same Efimov state 
(connected by the dashed line in Fig. [2| . The dashed line 
in Fig. [3] differs from the solid one by the scaling factor 
exp(7r/so) and shows the ratio al"'''"'^''/|a^''|. Note that 



FIG. 3: (Color online) Solid line: a, /\a_ \ vs 5, where n 
is the index of the Efimov state [see Eq. ([TJ]. Dashed line: 

(n + l) /I in) \ I / \ in) i, in), 

al /\aL \ = exp(7r/so)a, /\a- \- 



FIG. 5: (Color online) The coefficient C for the different 
three-body recombination rates in dependence of the mass 
ratio 5 — mxjm^. The inset shows C for larger values of 5. 




(a) 



FIG. 4: Diagrammatic representation of (a) the three-body 
recombination amplitude and (b) the elastic three-body scat- 
tering. Line patterns are the same as in Fig. [l] 



the scaling factor rapidly increases with 6 for (5 > 1 and 
one can conclude that a sequence of Efimov resonances 
is more likely to be seen in systems with smaller mass 
ratios. 



B. Three-body recombination for a > 

Let us now discuss the shapes of the inelastic loss res- 
onances and calculate the three-body rate constants in 
a heteronuclear system. We first consider the case of 
positive scattering length, a > 0, where the atoms can 
recombine into the shallow dimer and into deep dimers. 
The recombination into the shallow dimer can be related 
to the T-matrix element shown in Fig. |4];a). The event 
rate constant for inelastic scattering a is defined by the 
rate equation 



— n2 = 2— ni = -2anin2 , 
at at 



(10) 



where denotes the atomic number densities of the cor- 
responding species. 

The rate constant ag for recombination into the shal- 
low dimer is given by 



4m 



AD 



(11) 



with the dimer breakup momentum kjj — \J pi ad I p^a^^ . 
If deep dimers are present, their effect on the recombina- 
tion into the shallow dimer can be incorporated by an- 
alytically continuing the three-body parameter into the 
complex plane [18]. We thus make the substitution 



A — > Acxp(i?7*/so) 



(12) 



in Eq. Q, where 77* accounts for the effect of the deep 
dimers. nonzero value of 77 also generates the width of 
the Efimov trimers. By evaluating Eq. (11) numerically. 



we find that the known analytical formula for the three- 
boson case f24 simply acquires a new mass-dependent 
overall coefficient. The modified analytical formula is 
hence 

D(sin^[soln(a/ao*)] +sinh^77*) ^a'' 

OLs — C(0) 2 Al'jj 

sinh (ttsq + 77*) -f cos2[so ln(a/ao*)] mi 

where D = 1287r2(47r - 3\/3) and the mass-dependent 
coefficient is denoted by C((5). The parameter ag* gives 
the position of the minimum in the three-body recom- 
bination. The coefficient C(S) is shown in Fig. [5] The 

error in the extraction of C(5) from fitting Eq. (|13[) to 
our numerical results for is of order 10~^ for b< 2. 
For larger values of 8 the numerical extraction of C be- 
comes difficult because of a very large value of the scaling 
factor. To depict C{S) for i 5 > 2 we use the analytical 



formula (56) derived in Sec. IV 



Although our calculations in this section are conducted 
by varying the complex three-body parameter A, we 
present the results in terms of the practically relevant 
length parameters a_ < 0, a* > 0, and oq* > 0, and 
the dimensionless elasticity parameter 77*. The univer- 
sal theory predicts that the ratios a*/|a_| and |a_|/ao* 
depend only on the mass ratio (5. The former is shown 
in Fig. |3] and the latter is |a_|/ao» = exp(7r/2so) as is 
derived in Sec. IV This fixes the relative positions of all 
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the three-body loss features on both sides of the Feshbach 
resonance. 

The total rate of three-body recombination into all 
dimers (shallow and deep) for a > can be obtained 
from the optical theorem. It relates the imaginary part 
of the forward T-matrix element [shown in Fig. |4]jb)] for 
vanishing momenta to the event rate constant of inelastic 
scattering, a. This leads to the total recombination rate 
constant 

as + ad^ ImTi22^i22 = ST:aHmA{0, 0; 0) , (14) 

where A denotes the appropriately infrared subtracted 
amplitude 23]: 



Aip, k; E) = A{p, k; E) 
47r(l-H(5)2 



47r(l + (5) 



miap 
arcsin [1/(1 + 5) 



8a 
nil 



mip 

(1 -t- (5)2 arcsin [1/(1 -|- S)] - a/'J(2 + 6) 



Inp . 

(15) 



By subtracting Eq. (13) from Eq. (14) we find the rate 



constant for the recombination into deep dimers: 

D coth(7rso) cosh(?7*) sinh(77*) ha^ 



ad = C{S) 



sinh2(7rso + V*) + cos2[so ln(a/ao*)] mi 



,(16) 



where the coefficients C{S) and D are the same as in 
Eq. 

When So is not too small such that exp(27rso) ^ 1 (see 
Table m, the denominators in Eqs. ( 13 ) and ( 16 ) are prac- 
tically independent of a. In this case the a dependence 
of as and ad is simplified, and the corresponding expres- 
sions are known in the case of three identical bosons (see, 
for example, Ref. [1]). 



C. Atom-dimer scattering 

On the positive side of the Feshbach resonance (a > 
0) it is also possible to prepare an ultracold mixture 
of atoms and weakly bound dimers (see, for example, 
Refs. lElIlH])- An important observable in this case is 
the atom-dimer scattering length. Within our theory, it 
is given by 



a-AD = ~— — A 1 0, 0; - ^i^— ^ 



(17) 



27r V ' ' 2^a2 
and its universal dependence on a is parametrized by 

(Ci{6) + C2{S) cot[so ln(a/a,)])a , (18) 



a-AD 



where the coefficients Ci{S) and 6*2(5), calculated nu- 
merically, are shown in Fig. |6] Here we estimate the 
numerical error in the determination of Ci{S) and C2{S) 
to be of order 10"'^. 




FIG. 6: (Color online) The parameters Ci{S) and 6*2(5) in 
the expression for the atom-dimer scattering length, Eq. (fTsl . 



Efremov and collaborators have recently derived 
Eq. ( 18 1 for the atom-dimer scattering length in the 



Born-Oppenheimer approximation valid in the limit S — >■ 
[31]. For S = 0.081, corresponding to the ^Li-^'^Rb-^'^Rb 
system, our values for the coefficients C*i and C2 agree 
with the ones given in Ref. [31] to within 2-3% (see Ta- 
ble |l]) . However, we observe a stronger discrepancy in be- 
tween our value, sq = 1.523, and the Born-Oppenheimer 
result. So = 1.322, for this system [ST] . 

The effect of deep dimers on the atom-dimer scat- 
tering process can be incorporated by replacing a* — >■ 

exp(— i77*/so), equivalent to Eq. (12 1. At the scatter- 



ing threshold, the atom-dimer relaxation rate constant 
/3, defined by the rate equation 



d d 

-jinA = -rno = -prLATLD ■ 
at at 



(19) 



is given by [8] 

I3{E = -Bd) = -{inh/^iAD) ImaAD 

+ 1 sin"' [so ln(a/a*)] + sinh 77* mi 

(20) 

Furthermore, we can calculate the atom-dimer relax- 
ation rate constant above threshold. It is related to the 
inelastic atom-dimer scattering cross section by 



(3iE)^ 



(inol.) 
^AD 



(21) 



where k = ^2^iAD{E + Bd). The energy dependent in- 
elastic cross section is given by the difference of the total 
and the elastic cross sections, 
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FIG. 7: (Color online) The dimer relaxation rate constant /3 
in units of ha/nii for =0.1 and S = 0.471 as function of 
a/a*. The solid, short-dashed, long-dashed, and dot-dashed 
lines show /3 for E/Bd = —1, —0.95, —0.5, and 0, respec- 
tively. The double-dot-dashed line indicates the trajectory of 
the resonance maximum as the energy is increased from —Bd 
to zero. 



We can use this formula up to the dimer breakup thresh- 
old at k = kjj and thus map out the trajectory of the res- 
onance peak. It moves from a* at the scattering thresh- 
old, E — —Bd, to |a_| at the dimer breakup threshold, 
E — 0. For S < 3.475 the resonance peak moves to val- 
ues a > a*. Starting at S — 3.475, where we have exactly 
a*/|a_| — 1, it reverses this behavior and moves to val- 
ues a < a*. The peak height diminishes considerably 
with the energy. This effect is very large, especially for 
small values of S. For example, for 77* =0.1 and <5 = 0.1, 
the peak at = is smaller by a factor of 706 than the 
peak at E = —Bd- For (5 = 10 this ratio still is 19.1. At 
= 0, we find excellent agreement with the analytical 
formula 



P{E = -0) = n[S{5 + 2)f'^/{5 + If 

sinh(27rso) sinh(2?7^) ha 

X 7^ 

sinh (ttsq + 77*) + cos2[so hi(a/ao*)] "mi 



(23) 



which is derived in Sec. IV The peak position of (3 {E 



— 0) coincides exactly with the position of the maximum 
of the three-body recombination rate at threshold. 

In Fig. [7] we show numerical results for /3 for 77* =0.1 
and 6 = 0.471 corresponding to the K-Rb-Rb system ob- 
served in the Florence experiment [15] . The solid, short- 
dashed, long-dashed, and dot-dashed lines show /3 for 
E/Bd = —1, —0.95, —0.5, and 0, respectively. As the 
energy is increased toward the breakup threshold, the 
resonance height decreases strongly and the resonance 
becomes less pronounced. The double-dot-dashed line 
shows the trajectory of the resonance maximum as the 
energy is increased from —Bd to zero. As the energy is in- 
creased, the resonance position is shifted from a* toward 
larger values of a until it reaches its maximum value of 
2.04a* for E/Bd ~ —0.25. For larger energies, the res- 



onance position moves back to smaller values of a and 
reaches |a_| = 1.89a* at the dimer breakup threshold. 



D. Three-body recombination for a < 

On the negative side of the Feshbach resonance, shal- 
low dimers are absent and atoms can only recombine into 
deep dimers. The corresponding rate constant is again 
determined by using the optical theorem 



ad = ImTi22^i22 - STra^ W(0, 0; 0) 



(24) 



We have performed numerical calculations of ad for mass 
ratios S < 2, where the numerical accuracy is better than 
0.1%. Our results agree with the formula 



ad 



C{S) i:»coth(7rso)sinh(2r7*) ha'^ 
2 sin^[so ln(a/a_)] -f sinh^(?7*) TOi 



(25) 



where the coefficients C{S) and D are the same as in 
Eqs. (13) and (16). Equation (25) is derived in the next 



section. 



IV. ANALYTICAL APPROACH 

Relying on an analysis of the recombination process in 
configuration space, Macek et al. [32l |33], Petrov [34] . 
and Gogolin with co-workers [35' have obtained a com- 
pletely analytic expression for the three-body recombina- 
tion rate constant for identical bosons with a large pos- 
itive scattering length a when there are no deep dimers 
(i.e., r/* = 0). Braaten and Hammer generalized this 
result to finite ry* by making the analytic continuation 
K* — )■ K* exp(i77*/so) in the amplitude of this process [T]. 

In this section we would like to present an approach 
that allows one to calculate almost all zero-energy three- 
body observables analytically [in particular, a^, ad{a > 
0), and ad{a < 0)] for heteronuclear Bose-Fermi or Bose- 
Bose mixtures. We generalize the method introduced in 
Ref. |36j for fermionic mixtures and present it in more 
detail, also because it is conceptually different from the 
approaches of Refs. [5^1551155] . 

Let us consider the process of three-body recombina- 
tion into the shallow heteronuclear dimer. When the to- 
tal energy of the system, E, becomes smaller than the 
binding energy of the dimer, the three-body recombina- 
tion rate tends to a finite zero-energy value. In order to 
calculate this quantity we work in the coordinate repre- 
sentation. It is instructive to derive the real-space ana- 
log of Eq. ([4]) directly from the Schrodinger equation. At 
E — the latter reads 



E 

1=1,2 



2m2 



;7(|r-Ri|) 



2mi 



* = 0, 



(26) 
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where Ri are the coordinates of the bosons and r is the 
coordinate of the third atom. For the interspecies inter- 
action we use the zero-range Fermi pseudopotential 



C/(|r-Ri|)- 



-S{r-Ri) 



d 



d\r-Ri\ 



|r-Ri|-)- 



(27) 

Let us temporarily adopt the units h — 2fj, ^ 
a — 1, separate the center-of-mass motion, and in- 
troduce the new coordinates x = [Ri sin (j) — R2 + 
r(l — sin (/>)]/ cos </) and y = r — Ri, where is de- 
fined in Eq. ([8|. The bosonic symmetry condition, 
\I'(Ri,R2,r) — \I'(R2, Ri, r), in the new coordinates 
reads 

^'(x, y) = X sin (f) + y cos (j), x cos (j) + y sin (j)) , (28) 



and Eq. ( 26 ) takes the form 



(-V2-V^)vI/(x,y) = /(x)5(y) 

+/(— X sin (j) + y cos (/))(5(x cos - 



ysin0), (29) 



where we define 



/(x) 



-An \iind[y^{x,y)]/dy. 



From Eqs. (29) and (30) it is easy to see that 



*(x,y) 



>/(x)(l/j/-l)/47r. 



(30) 



(31) 



Thus, the function / describes the motion of an atom 
relative to the other two atoms when they are on top of 
each other. In some cases it is useful to consider / as an 
atom-dimer wave function. Indeed, for large x expres- 
sions 5* oc 4>b{y) e'yi'p{ix) / X and f{x) oc (Stt)-'^/^ exp(ia;)/a: 
equivalently describe the relative outgoing motion of the 
products of three-body recombinatio n. The coefficient 
(Stt)^/^ is obtained by comparing Eq. (31 ) and the small- 



y asymptote of the normalized molecular wave function 



0b(y) = exp(-2/)/V2^y — > [l/y - 1)/V27r. (32) 
■y-i-o 

Therefore, in order to calculate the three-body recombi- 
nation rate constant we first derive and solve an equa- 
tion for /, then separate the large- a; asymptote / oc 
exp(z2;)/a;, and finally relate the coefficient in front of 
it to the three-body recombination rate constant. 

Using the Green's function G{X) = l/An^X^ of 



Eq. ( 29 ) we can express ^ through /; 



*(x,y) = *o + y [G(V(x-x')2 + y2) 
-f'G(v/(x + x' sin(/))2 + (y-x' cos <?!))2)]/(x')d^a:', (33) 



where ^0 is a solution of the homogeneous Euler equation 
(— Vx — Vy)5'o = without singularities. In the case of 
three-body recombination, 4'o is a correctly normalized 
wave function of three free atoms. We consider the atoms 



in a volume V in a state where the two bosons are not 
in the same quantum state (cold thermal gas). Then, in 
the region relevant for the recombination (x < 1, y < 1), 
we get 



'2V' 



(34) 



We point out that vj/g should be set to zero for the prob- 
lem of atom-dimer scattering just below the breakup 
threshold {E = —0) as in this case the atoms cannot 
move freely at large distances. 



The equation for / is obtained by substituting Eq. ( 33 ) 
into Eq. (30). We write it in the form 



{L - l)/(x) - 4^^'o, (35) 
where the integral operator L is defined as 

L/(x)=4^|{G(|x-x'|)[/(x)-/(x')] 

-G{y/x'^ + a;'2 + 2xx' sin (t>)f{'x!)}dPx'. (36) 

The operator L conserves angular momentum, and 
Eq. ( 35 ) can be written as a set of uncoupled equations 
for each spherical harmonic of /(x). For the three-body 
recombination problem at hand we look for a spherically 
symmetric solution /(x) = f{x). Then, Eq. (36) reduces 
to 



(a;2 — x'2)2 



{x^ + x'^y - Ax 



2x'2 sin 



x'^dx'. (37) 



The integrals in Eqs. (36) and (37) are taken as principal 
values [37]. 

Let us now discuss the structure of possible solutions of 
Eq. ( 35 ) . Obviously, f{x) is a sum of a particular solution 
of the inhomogeneous Eq. ( 35 ) and a general solution of 
the homogeneous equation 



(i-l)x(x) = 0. 



(38) 



Physically, Eq. (|38j) describes the atom-dimer channel 
just below the dimer breakup threshold (^q — 0). There- 
fore, at distances a; 3> 1, the function x(a;) is a linear 
combination of exp{ix)/x and exp{—ix)/x. 

In order to understand the short-distance behavior of 
x(a:), note the following property of the operator L: 



Lx" = X{v)x''- 



(39) 



for any complex exponent v in the region —3 < Re(i^) < 
1, which is, in fact, the region of convergence of the inte- 
gral on the left-hand side of Eq. ( 39 ) 38J . The function 
\{v) is given by 



2 sinfi 



1) 



sin(2(/)) cos(7ri'/2) 



(40) 
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and has two complex conjugate roots, = — 1 ± isq, 
where Sq is a real number satisfying Eq. ([7|. At short 
distances the operator L in Eq. (38) dominates over 1, 
and any solution of this equation should be a linear su- 
perposition of X oc x~^~^^^'' and its complex conjugate. 
From now on we use the notation x for the solution of 
Eq. (38) with the following asymptotes: 



Ax"- 



^ — l^ix+icr — h 



a; < 1 
X > 1, 



(41) 



The rule of going around the pole in the second inte- 
gral and the numerical coefficient in front of the second 
term in the square brackets on the right-hand side of 
Eq. (46) regulate the entry of xi^): which can be ar- 
bitrary, into the particular solution go{x). Using this 
freedom, we choose these parameters in such a way that 
go{x) does not contain oscillating terms proportional to 
^.-i-l-iso small X. Direct calculation shows that in 
the limit x — > the right-hand side of Eq. ( 46 ) equals 
.9o(a^) ~ [/o x(^)dz]^/27ra; to the leading order in x. 
On the other hand, according to Eq. (39 1, the same 



^^^'''r^J^l':'''^'^}''™^^^^^^^ quantity in the same limit can be writtm as go{x) = 

^ = 1/[A(— l)a;], which leads to 



bers. The physical solution of Eq. (38) (i.e., the one cor- 
responding to a given three-body parameter) is expressed 
as 



Xe{x)^e^\{x) + e~^'x*ix), 



(42) 



where 6 is the three-body parameter (a complex number 
with imaginary part 77*). 



The normalization in Eq. (41 1 is chosen such that 



{px{px)\p'x{p'x)) = / px{px)p'x{p'x)x dx = 2tt6{p-p'). 
Jo 

(43) 



The first equality in Eq. ( 43 1 is our definition of the scalar 



product (note the absence of the complex conjugation), 
and the second equality follows from the fact that px{px) 
and p'xip'x) are eigenfunctions of the symmetric opera- 
tor L corresponding to the eigenvalues p and p' . They 
are orthogonal for p ^ p' and their scalar product in the 
vicinity oi p = p' can be worked out in the same way as 
in Ref. [35] (see Sec. 21 there). A simple change of the 
integration variable in Eq. ( 43 ) leads to the completeness 
condition 



p^x{px)x{px')dp — 2nS{x — x')/x'^ 



(44) 



Equations ( 43 ) and ( 44 ) allow us to construct the inte- 
gral operator (i — 1)^^ needed to solve Eq. (35 1. In order 



to avoid problems with divergence of the corresponding 
integrals let us introduce an auxiliary function ga{x) re- 
lated to f{x) by 

fix) = 47r^o[-l - m/x + A(0)A(-l).go(a;)]. (45) 



Subs titu ting this expression into Eq. ( 35 ) and using 



Eq. (39), we find that go{x) satisfies the equation {L 



l)go{x) = X . Applying the operator (L — 1) 
we obtain the following particular solution: 



to X 



1 



9o(x) 



X{z)dz 



lixx Jo 

°° _x{y)ydy 

Uo 



2'!Tixx{x) 



y — X — iO 1 — exp(— 27rso) 
where the first integral is defined as 



xiz)dz = lim 



+0 



X{z)z''dz. 



(47) 



{L-iy^x- 

the result 



L- 



X{z)dz - V2^/A(-l). 



(48) 



Another consequence of our choice of the particular so- 
lution ( 46 ) is that removing the oscillating terms from 



5o(x) makes it real, since any imaginary part of go would 
necessarily be a solution of the homogeneous Eq. ( 38 1 . 



Therefore, g^ would have oscillations at short x, the ab- 
sence of which we have ensured. Clearly, the function 



/ obtained by virtue of Eq. (45) is also real. Moreover, 



property (39 1 ensures that / = o(l) at small x (i.e., its 
Taylor expansion starts with at least). Therefore, this 
solution of Eq. (35) is not sensitive to the short-range 



physics and does not depend on the three-body parame- 
ter. 



Integrating Eq. ( 46 ) in the limit x 3> 1 we get 



50 (x) 



2ti 



i{h + vrso)] 



>oo sinh(7rso) V A(— 1) x 

(49) 

It can be real only ii h — — ttsq [note that A(— 1) < 0]; 
compare Ref. [52] . 

Finally, the result that we are interested in is the linear 
combination 



feix) = fix) +-fxeix), 



(50) 



where the complex number 7 is chosen such that feix) 
contains only an outgoing wave at large x (corresponding 
to an atom and a dimer flying apart after the three-body 
recombination event). This condition gives 



7 



7r*oA(0)^27rA(-l) 
sinh(7rso) cosh(7rso — iO) 



(51) 



Keeping only the relevant oscillating term at large x, we 
obtain 



fe{x) 



-> iA'y sin 9 sinh(7rso) exp(ix -|- ia) jx. (52) 



(46) Equation (l52| together with Eqs. (pll) and (pSl) gives 



us the atom-dimer outgoing flux. Indeed, the large-x 
asymptote / = ^exp(zx)/x, where ^ is any complex am- 
plitude, is accompanied by the flux |^P$oo7 where 



2 X (87r)~i X (47r) x 2 



(53) 
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Here we have explicitly written out the following factors: 
the factor of 2 reflects the two symmetric possibilities 
of forming the dimer (corresponding to the interchange 
Ri ^ R2), the factor of (Stt)"^ arises from the rela- 
tion in between 5* and / [see the discussion preceding 
Eq. (32)], the factor of Att is the solid angle in the out- 
going atom-dimer channel, and the last factor of 2 is the 
atom-dimer relative velocity in the x, y coordinates. The 
three-body recombination rate constant is obtained 
by taking the squared modulus of the prefactor in front 
of exp{ix + ia)/x in Eq. ( |52| ) and by multiplying it by 
$00 5 by the factor of 1/2, reflecting the fact that the 
number of pairs in the gas of species 2 is n|/2, and by 
the factor ha'^/2fi in order to restore the original physical 
units. We should also mention that the nine-dimensional 
volume is taken to be a unit volume in the original sys- 
tem of coordinates {r, Ri,R2}. In the new coordinates 
{x, y, Rem.}, where Rem. is the center-of-mass coordi- 
nate, this volume is — cos~^ 4>. The final result for 
the three-body recombination rate constant reads 

327r3cos3<^A2(0)|A(-l)| 
sm (h 



sin^[so ln(a/a*o)] + sinh^ 77* ha^ 
sinh^(7rso + V*) + cos2[so ln(a/ao*)] nii 



(54) 



where we have expressed the three-body parameter 9 
through the original physical units: 



So ln(a/a*o) + iv* 



(55) 



Formula ( 54 1 is in excellent agreement with our numerical 
results. Comparison with Eq. ( 13 ) leads to 0D] 



C{S) 



(1 + (5)2 arcsin [1/(1 + S)] - ^S{2 + S) 



(56) 



2(47r - 3^3) 

As explained in Sec. |III| the constant of the three-body 
recombination into deep dimers can be derived from the 
optical theorem. The result for Ud is given by Eq. ( |16| . 
Here we would like to show how one can derive this result 
by using the method of this section. In contrast to the 
recombination into shallow states we now have to look at 
the balance of the incoming and outgoing fluxes of atoms 
corresponding to the short-distance asymptote of fe{x) 
given by Eq. ([50| 

fe{x) > jxe{x) = j{Ae''x-^+''° +A*e'''x''-''^'), 

(57) 

In analogy with $007 let $0 denote the number of atom 
triples disappearing at the origin (x = 0, y = 0), pro- 
vided the function / takes the form of the incoming wave 
^-i-isa TffiiYy ynji; weight. With this definition, the re- 



combination rate constant follows from Eq. (57) 
ad = {hay2^iMHo\A\^smh{2r^,). 



(58) 



The product <I>o|A|2 can easily be found from definition 
( 41 ) by equating the fluxes at a; — > and at x — >■ 00 and 



using Eq. (53): 



$o|^|^ = 2<I>oo sinh(27rso) = 4sinh(27rso) 



(59) 



Substituting Eqs. (59) and (51) into Eq. (58), we obtain 



Eq. ( 16 ) exactly. We point out that it is in principle pos- 
sible to obtain $oi and, therefore, by substituting 
the expression / = x~^~^''° into Eq. (33) and calculating 



the incoming flux from the resulting wave function \l/ in 
the six-dimensional conflgurational {x, y}-space. 

Our analytical approach can also be used to derive 
on the negative side of the resonance (a < 0). In this 
case we use the units h — \a\ — 2fi = 1, and some equa- 
tions described above should be modifled accordingly. In 
particular, Eq. (35) reads 



(L-f l)/(x) =4^«'o 



(60) 



and we now introduce an auxiliary function ^0 related to 
/by 



fix) - 4^*o[l - A(0)/a; + X{0)X{-l)g„{x)] 



(61) 



where go satisfles {L + l)go{x) ^ x We write the 
solution in the form 

goi^) ^ ^ r xi^)dz r (62) 

Sttx Jo Jo V + x 

and integrating it in the small-x limit we get the asymp- 
tote 



fix) 



27r*oA(0)v/27rA(-l) 
sinh(7rso) 



Ax-^+''°. 



(63) 



The function / is a solution of Eq. (60), but its oscil 



lations at small x do not have (in general) the correct 
phase imposed by Eq. (41 ). This difficulty is resolved by 



observing that /* also satisfies Eq. (60). The correctly 
behaving solution reads 



fe{x) 



exp{i0)f{x)+cxpi-te)r{x) 
exp{i9) + exp{—i9) 



(64) 



and by subtracting the outgoing flux from the incoming 
one at small x we obtain the result, cf. Eq. (|25|), 



ad(a < 0) = 167r3cos^(/)A2(O)|A(-l)|coth(7rso)/sin0 
sinh(277*) ha^ 



cos2[so ln(|a|/ao*)] + sinh^ 77* mi 



Equation ( 65 ) also gives the ratio 

|a_|/ao* = exp(7r/2so). 



(65) 



(66) 



In other words, the maxima of ag and ad (a > 0) and the 
maxima of ad (a < 0) are symmetric with respect to the 
center of the Feshbach resonance. 

Let us now return to the case a > and discuss some 
properties of atom-dimer collisions just below the dimer 
breakup threshold (E = —0). Namely, by substituting 
the large- a; asymptote of x into Eq. ( 42 ) one readily ob- 



tains the s-wave contribution to the scattering S'-matrix 



5*0 



2ia 



cosh(7rso + i^) / cosh(7rso — i6 



(67) 
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From this expression, one can calculate the atom-dimer 
inelastic rate constant at zero total energy: 



f3{E = -0) = vr cos^(/i cot (/i 

sinh(27rso) sinh(277*) 



ha 



sinh (ttso + r?*) + cos^[so ln(a/ao*)] 'm-i 



(68) 



see also Eq. (23). Remarkably, (3 {E ~ —0) also reaches 
its maximum at a = |a_|, and for small sq the peak can 
be quite narrow. We conjecture that this behavior is due 
to the Efimov state that crosses the atom-dimer thresh- 
old at a — and then exists as a scattering resonance. 
Our numerical results support this conjecture. We find 
that for mass ratios S < 3.475 the resonance peak moves 
to values a > a*, while it moves to values a < for 
6 > 3.475. The resonance peak then hits the three-atom 
threshold at a = |a_| exp(— tt/sq)- In Fig. [tI we have 
shown /? for the specific case 77* — 0.1 and = 0.471 
corresponding to the K-Rb-Rb system observed in the 
Florence experiment [T5]. The trajectory of the scatter- 
ing resonance is given by the double-dot-dashed line in 
the figure. For smaller sq leading to a larger scaling fac- 
tor, we conclude that the center of this Efimov scattering 
resonance can travel quite far from the value a = a.^ at 
E = —Bd to the value a = |a_| exp(— tt/sq) at E — 0. We 
have verified this behavior numerically. It is thus worth 
investigating this scattering resonance in the vicinity of 
the atom-dimer threshold experimentally. 



V. COMPARISON TO EXPERIMENT 

There are several experiments on heteronuclear Bose- 
Bose and Bose-Fermi mixtures, to which our results are 
directly applicable (in the universal limit). In Table |lj 
we present the universal predictions for some combina- 
tions of alkali isotopes being investigated at the moment 
and interesting from the viewpoint of Efimov few-body 
physics. We sort them by the value of the scaling factor. 



A. The ^"K-'^'Rb mixture 

Zirbel et al. at JILA recently studied weakly bound 
fermionic ''"K-^^Rb molecules and their stability in colli- 
sions with atoms near a wide (open-channel-dominated) 
heteronuclear Feshbach resonance at Bq = 546.7 G [55]. 
In particular, they measured the atom-dimer relaxation 
rate for collisions of these dimers with Rb atoms as a 
function of a. The corresponding data (see Fig. 2 in 
Rcf. [1^) can be fit very well with our Eq. (20), where 
the fitting parameters are a* — 200 ± 50 asohr and 
r|^. — 0.05 ± 0.02. In the same work, the authors also 
measured the three-body recombination rate constant on 
both sides of this Feshbach resonance (i.e., in the same 
universal region) . We fit their results with Eq. ( 25 ) on 



the negative side of the resonance and with the sum 
Us + ad given by Eqs. ( 13 1 and ( 16 ) on the positive 
side [H]. A good agreement is achieved if we choose 
a* = 300 ± 100 QBohr (for this mass ratio a_ = —1.96a*, 
see Table |l| and the same 77* as above. These parameter 
values lead to a peak of the three-body recombination at 
a = a_ ~ — 600 ttBohr- Although in Ref. [2S] the peak 
has not been identified, the overall shapes of /3 and a 
measured for this particular Feshbach resonance indicate 
that it is worth performing a more detailed measurement 
of the three-body loss rate around this value of a. 



B. The **^Rb-^'K mixture 

The group of Inguscio and Minardi in Florence inves- 
tigated a Bose-Bose mixture of ^^Rb and ''^K ^5j. They 
observed three loss resonances by scanning the scatter- 
ing length and monitoring the population dynamics of 
the species in the vicinity of each of the resonances. For 
negative scattering length, they identified a K-Rb-Rb res- 
onance at a = —246 a Bohr and a K-K-Rb resonance at 
a — —22000 aBohr- The third resonance is observed at 
the positive scattering length a = 667 aBohr and at- 
tributed to enhanced atom-dimer scattering in the K- 
Rb-Rb three-body system. This process is assumed to 
contribute to three-body losses through multiple rescat- 
tering processes (see also Ref. tl2j). An independent con- 
firmation of this resonance in a system prepared directly 
out of K-Rb dimers and Rb atoms would be desirable. 
Unfortunately, in contrast to the JILA experiment, the 
dimers are bosonic and their short lifetime can make such 
a confirmation difficult [42j . 

The interspecies van der Waals length in the K-Rb 
system is ivdw = '^'^a.Bohr, such that these resonances 
should be within the range of validity of the universal 
theory. Assuming that the observed K-Rb-Rb features 
are due to Efimov resonances, one can extract the ra- 
tio a,t/|a_| — 2.7 from the Florence experiment, whereas 
our theory predicts a,t/|a_| = 0.52. The discrepancy 
can be attributed to the effective range corrections. In 
particular, one should be careful with the feature at 
a = —246 ttBohr, which is not too large compared to 
the van der Waals length. Besides, if we believe in the 
"rescattering" nature of the positive-a resonance, one 
should take into account a finite-energy shift of the po- 
sition of the atom-dimer scattering resonance. Indeed, 
even at zero temperature, dimers formed by three-body 
recombination collide with stationary Rb atoms at the 
finite collision energy [m'j^ij/{mK + '2mjib)'^]Bd ~ O.ldBd- 
In Fig. [Tj we have shown numerical results for /3 for 
energies from the scattering threshold, E = —Bd, up 
to the breakup threshold, E — 0, using — 0.1 and 
6 = 0.471 corresponding to the K-Rb-Rb resonance at 
a ~ 667 aBohr- At E = 0, the resonance peak is only 7% 
higher than the minimum value of /3 and the resonance is 
almost completely washed out. Moreover, the peak value 
of /3 at E' = is a factor of 300 smaller than at E = —Bd- 
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■^Li-Cs-Cs 


^Li-Rb-Rb 


■'Li-Rb-Rb 


^"K-Rb-Rb 


"K-Rb-Rb 




5 


0.053 


0.069 


0.081 


0.460 


0.471 


2.21 


So 


1.850 


1.635 


1.523 


0.6536 


0.6444 


0.2462 


exp(7r/so) 


5.465 


6.835 


7.864 


122.7 


131.0 


348000 


a,/\a-\ 


0.13 


0.16 


0.18 


0.51 


0.52 


0.91 


C{5) 


0.072 


0.068 


0.066 


0.037 


0.037 


0.015 


Ci{5) 


2.54 


2.33 


2.22 


1.14 


1.13 


0.94 


C2{5) 


2.52 


2.5 


2.47 


2.08 


2.07 


1.30 



TABLE I: Universal parameters for various heteronuclear mixtures. The isotopes of rubidium and cesium are «^Rb and "^'Cs. 
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FIG. 8: (Color online) aa vs a for the "^K-^^Rb-^^'^Rb system 
assuming a_ = —246 aBohr and three values of 77,: r;* = 0.12 
(dashed line), 77, = 0.01 (solid line), and 77* = 0.4 (dash- 
dotted line). Data points indicated by diamonds are taken 
from Refs. [I5l|42]; see text. 



Of course, the explicit numbers depend on the value of 
77*, but one should not exclude the possibility of other 
explanations of the positive-a feature. 

Aside from detecting the positions of the resonances, 
it is desirable to compare the actual shapes of the a- 
dependence of the loss rate with the theoretical predic- 
tions, especially on the positive side of the resonance, 
where as -f ad is rather smooth. So far, the three-body 
recombination rate in the ^^Rb-'^^K mixture has been 
measured for two values of a < [THlllS], and the compar- 
ison with our calculation is rather inconclusive. Figure [8] 
shows Ud for the K-Rb-Rb resonance at a = —246 asohr 
where only the recombination into deep dimers can occur. 
The dashed curve is calculated using the value 77* = 0.12 
suggested in Ref. [15J. The data point close to the reso- 
nance is taken from Ref. [13 US] , whereas the one farther 
away from the resonance gives an upper limit of the re- 
combination rate [32]. In order to illustrate the sensitiv- 
ity of the result to 77*, we also show it for 77* = 0.01 (solid 
line) and 77* = 0.4 (dash-dotted hue). The discrepancy 
between the measured recombination rate at the reso- 
nance and our result for r]^ = 0.12 is about one order of 
magnitude. In order to understand its origin, more mea- 



surements around the resonance position are required. 
Such data would allow for a more precise determination 
of the width parameter 77* and of the resonance shape 
predicted by the universal theory. 



C. Future experiments 

The Tiibingen group of Zimmermann recently stud- 
ied the ■^Li-^'^Rb and ^Li-^'^Rb mixtures. They identified 
and quantified several interspecies Feshbach resonances 
in both of them [331 135] and have reached quantum de- 
generacy |46j . These mixtures are characterized by quite 
small mass ratios and, therefore, small scaling factors, 
which is favorable for observing the discrete scaling in- 
variance. Another very good candidate for studying the 
Efimov effect with even smaller scaling factors is a mix- 
ture of ^^^Cs with either isotope of lithium (the ^Li-^^^Cs 
mixture was created in Heidelberg [AT ) . The universal 
parameters for these mixtures can be found in Table [l] 
Predictions for the three-body recombination and atom- 
dimer relaxation rates can be obtained from Eqs. ( |13[ ), 
diel), pOl, (1231), and 



VI. SUMMARY AND CONCLUSION 

In this paper, we have calculated the three-body loss 
rates in heteronuclear mixtures of atoms for the case of 
large scattering length between the unlike atoms. We 
have analyzed this problem using two complementary 
methods. 

First, we have formulated a universal effective field the- 
ory for this system and derived momentum-space integral 
equations for the trimer energies and the atom-dimer 
scattering amplitude. From an analysis of the bound- 
state equation we have calculated the ratio of the reso- 
nance positions a,t/|a_| as a function of the mass ratio, 
5. Moreover, we have calculated the three-body recombi- 
nation and atom-dimer relaxation rates numerically. We 
have provided semianalytical expressions for the rate con- 
stants of three-body recombination into shallow and deep 
dimers as a function of the interspecies scattering length 
a and the Efimov width parameter 77*. Furthermore, we 
have calculated the atom-dimer relaxation constant from 
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the scattering threshold aX E = —Bd up to the dimer 
breakup threshold at i? = 0. 

Second, we have carried out an analysis of the recom- 
bination and relaxation process in configuration space. 
We have generalized the method developed in Ref. [SS] 
to heteronuclear bosonic mixtures and obtained analytic 
expressions for the recombination and relaxation rates at 
E = Q. We find excellent agreement of these expressions 
with our numerical results from the momentum-space in- 
tegral equations. 



The expressions in Eqs. (13), (16l, (25), and (56) fully 



determine the three-body recombination rates for het- 
eronuclear bosonic mixtures with resonant scattering be- 
tween the unlike atoms in the universal zero-range theory. 
The atom-dimer relaxation rates at _E = —B^ and E = Q 



are given by Eqs. (18), (20), and (23). These equations 



are universal and can be used to analyze experimental 
data for any combination of atoms with the range appli- 
cability of the universal theory. 

In Ref. [IH], D'Incao and Esry give a general functional 
dependence of the recombination rates on the scatter- 
ing length for all possible combinations of bosons and 
fermions. This includes the case of two identical bosons 
and a third atom with J = which we address here. We 
agree with their expressions for ad in the case a < and 
for p. For a^, our general form (54) does not agree with 
their result. The proportionality of to sin^(so lna+(/)3) 
where (t>^ is a short-range phase [HI emerges only if 
exp(27rso) 3> 1 and the expression ( |54[ ) can be simpli- 
fied (cf. Ref. Ij). This is the case for small mass ra- 
tios 5. Moreover, our prediction for the dependence of 
a., on 5 [see Eqs. (13) and (56)] differs from the result 



as cx [(5(2 + (5)]3/2aV(l + SfJ^i obtained in Ref. [48]. 

We have applied our results to some heteronuclear mix- 
tures in ongoing and planned experiments. We find good 
agreement between theory and the JILA experiment |28j 
that investigated ^°K-^^Rb molecules and their stabil- 
ity in collisions with atoms near a wide heteronuclear 
Feshbach resonance at i?o = 546.7 G. For the recent ex- 
periment by the Florence group which uses a mixture of 
''^K and ^'^Rb atoms [15, 43j, we observe moderate dis- 
crepancies between theory and experiment. We obtain 
a*/|a_| = 0.52 for the resonance positions while the ex- 



perimental ratio is a*/|a_| = 2.7. Because neither the 
effective range corrections nor the experimental errors of 
the ratio are known accurately, no definite conclusion can 
be drawn at the moment. In particular, our analysis of 
atom-dimer relaxation suggests that explanations should 
be considered other than an Efimov resonance for the 
feature at a = 667 asohr that was used to extract the 
value of a*. 

Using the value 77* = 0.12 extracted in Ref. [TS], we 
find that the calculated recombination rate at the reso- 
nance is about one order of magnitude too small. Using 
smaller values of 77*, the size of the experimental rate 
can be reproduced. In order to resolve this discrepancy, 
more measurements around the resonance position are re- 
quired. Currently, there are only two data points and 77,^ 
cannot be determined accurately. Additional data would 
allow for a more precise determination of rj^ and allow for 
a test of the resonance shape predicted by the universal 
theory. 

Finally, we have calculated the universal parameters 
determining the three-body loss rates for various other 
mixtures and have summarized them in Table ID Extend- 
ing earlier work by D'Incao and Esry [351 HH], our pre- 
dictions lay the theoretical ground for the experimental 
observation of Efimov physics in heteronuclear mixtures. 
They should be useful for planning and analyzing future 
experiments. 
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